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ABSTRACT 

Context. Recent observations of hard X-rays and very high energy gamma-rays from a number of young shell type 
supernova remnants indicate the importance of detailed quantitative studies of energy spectra of relativistic electrons 
formed via diffusive shock acceleration accompanied by intense nonthermal emission through synchrotron radiation and 
inverse Compton scattering. 

Aims. The aim of this work was derivation of exact asymptotic solutions of the kinetic equation which describes the 
energy distribution of shock-accelerated electrons for an arbitrary energy-dependence of the diffusion coefficient. 
Methods. The asymptotic solutions at low and very high energy domains coupled with numerical calculations in the 
intermediate energy range allow analytical presentations of energy spectra of electrons for the entire energy region. 
Results. Under the assumption that the energy losses of electrons are dominated by synchrotron cooling, we derived 
the exact asymptotic spectra of electrons without any restriction on the diffusion coefficient. We also obtained simple 
analytical approximations which describe, with accuracy better than ten percent, the energy spectra of nonthermal 
emission of shock-accelerated electrons due to the synchrotron radiation and inverse Compton scattering. 
Conclusions. The results can be applied for interpretation of X-ray and gamma-ray observations of shell type supernova 
remnants, as well as other nonthermal high energy source populations like microquasars and large scale synchrotron 
jets of active galactic nuclei. 
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ray emission. Nevertheless, the inverse Compton origin of 
gamma-rays remains an alternative option, provided that 
the magnetic field in the gamma-ray production region does 
not exceed 10/xG (Aharonian et al. I2006p . 

Power-law distributions of relativistic particles are read- 
ily formed at astrophysical shocks, in particular in young 
shell type SNRs, through the so-called diffusive shock ac- 
celeration (DSA) mechanism (Krymsky 119771 Bell 119781 
Axford et al. I1977[) . The maximum energy of accelerated 
particles and the shape of the spectrum around and be- 
yond the maximum energy (the so-called "cutoff region") 
is determined by the competition between the acceleration 
and escape rates, as well as, in the case of electrons, by ra- 
diative (synchrotron and inverse Compton) energy losses. 
In this regard, the most important information about the 
accelerator is contained in the energy distribution of parti- 
cles in the cutoff region. The spectrum of relativistic elec- 
trons in the cutoff region is formed under conditions when 
particles gain energy during the shock crossing and lose en- 
ergy simultaneously. In young SNRs, the radiative cooling 
of electrons results in formation of highest energy tails of 
synchrotron X-rays and inverse Compton gamma-rays, at 
> 1 keV and > 10 TeV, respectively. Thus, the compari- 
son of hard X-ray and TeV gamma-ray observations with 
the spectral features predicted by the DSA model in these 
energy bands gives direct information about the key param- 
eters characterizing the process of particle acceleration. In 



1. Introduction 

The nonthermal X-ray emission detected from a number 
of young shell-type SNRs (for a review see e.g. Vink 2006) 
is generally interpreted as synchrotron radiation of rela- 
tivistic electrons accelerated diffusively by shock waves to 
multi-TeV energies (see e.g. Drury et al. I2001|) . Two other 
(alternative) radiation mechanisms - the bremsstrahlung 
of subrelativistic electrons and the inverse Compton (IC) 
scattering of moderately relativistic electrons - are not suf- 
ficiently effective to explain the observed X-ray fluxes. The 
discovery of TeV gamma-rays from young SNRs, in par- 
ticular from Cas A ( Aharo nian et al. [2lKjT]), RX J1713.7- 
3946 (Enomoto et al. [20021 Aharon ian et al. 12004]) a nd RX 
J0852.0-4622 (Katagiri et al [203(751 Aharonian et al. 
provides unambiguous evidence of acceleration of particles 
(electrons and/or protons) to energies 100 TeV and beyond. 
Unfortunately even the very high quality morphological and 
spectrometric studies of young SNRS in TeV gamma-rays 
performed with the HESS array of imaging atmospheric 
Cherenkov telescopes do not allow robust conclusions con- 
cerning the origin of gamma-rays. For example, in the case 
of the best studied gamma-ray emitting SNR, RX J1713.7- 
3946, the hadronic model can explain satisfactorily both the 
overall energetics and the spectral features of TeV gamma- 



Send offprint requests to: V.N. Zirakashvili 



2 



V.N. Zirakashvili, and F.Aharonian: Energy spectra of electrons accelerated in supernova remnants 



1 








Ui 




upstream 

/ 


\ 

S 

X 


downstream 


/ 

/ 

/ 

/ 

/ 

/ 

/ 


N 

X 


. N(x) 






u 2 








X 





Fig. 1. The shock speeds in the upstream and downstream 
regions (solid lines) and the spatial distributions of elec- 
trons (dashed lines), accelerated at the plane shock. 



this regard, analytical presentations of the energy distri- 
butions of high energy electrons and the spectra of their 
synchrotron and IC radiation components provide effective 
tools for studies of nonthermal processes in SNRs. 

So far, analytical solutions for the electron spectra at 
the plane shock have been derived only for a special case of 
energy-independent diffusion coefficient (Bulanov & Dogiel 
[T9791 Webb et al. [T984l Heavens & Meisenheimer [T957|) . 
At the same time, for typical shock speeds of several thou- 
sand km/s in young SNRs, the particles, in particular elec- 
trons, have a chance to be boosted to multi-TeV ener- 
gies only if the acceleration proceeds in the most effective 
way, the so-called Bohm diffusion limit. This implies strong 
energy- dependence of the diffusion coefficient, D(E) oc E. 
Therefore, it is important to derive analytical solutions for 
the electron spectra for a more realistic case of diffusion. 

In this paper we present exact asymptotic solutions for 
the high-energy tails of distributions of electrons formed at 
the plane shock, assuming arbitrary energy-dependencies 
of both the diffusion coefficient and the energy lose rate of 
electrons. The spectrum of electrons in the low-energy do- 
main where the radiative cooling can be neglected, is well 
known from the solution of the transport equation char- 
acterizing the acceleration at the shock without losses. On 
the other hand, assuming that the cooling of electrons at 
very high energies is dominated by synchrotron losses, we 
obtained simple analytical expressions for the exact asymp- 
totic solutions in the high energy domain. In order to extend 
the analytic presentations to a broader energy range, we 
should "glue" the asymptotic solutions applicable at very 
low and very high energies, using numerical calculations 
performed at intermediate energies. We treat separately the 
electron distributions at the shock, as well as downstream 
and upstream of the shock, and calculate the spectra of 
synchrotron and IC radiation in these regions. We present 
these spectra in forms of simple analytical approximations 
which appear significantly different from ones often used 
in the literature for the fits of spectral measurements and 
subsequent interpretations of the X-ray and TeV gamma- 
ray observations of SNRs in terms of synchrotron and IC 
radiation of multi-TeV electrons. 



2. Asymptotic form of the electron spectra 

Below we consider acceleration of electrons by the plane 
shock. The upstream plasma moves with a speed u = u\ 
towards the shock from — oo of the x axis (see Fig. I). Then, 
the downstream speed is u 2 = wi/c, where a is the shock 
compression ratio. The electron momentum distribution 
N(p, x) obey the following equation written in the shock 
frame upstream and downstream of the shock: 



d dN dN 1 8 NAr „ 
— D— - u— + — —p 2 b{p)N = 
ox ox ox p z op 



(1) 



Here b{p) is the energy lose rate of particles (b > 0). 
Hereafter we assume that the synchrotron losses domi- 
nate over the escape of electrons from the system. The 
electron momentum distribution N(p) is normalized to 
n = An J p 2 dpN(p), where n is the number density of ac- 
celerated particles. 

Boundary condition at the shock front, x — 0, can be 
written as 



Ul 



P dN 
Is dp 
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dN 
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x=+0 



dN 
dx 



(2) 



x=-0 



where "f s = 3ct/(ct — I) is the power-law index of particles 
accelerated at the absence of energy losses, No(p) is the 
electron distribution at the shock front, D\ and D 2 are 
the diffusion coefficients upstream and downstream of the 
shock, respectively. 

Analytical solutions of Eqs.(I) and (2) can be obtained 
for simplified conditions, e.g. at low energies when the 
radiative energy losses can be neglected, or in the case 
of energy- independent diffusion-coefficient. On the other 
hand, it is possible to derive asymptotic solutions for an ar- 
bitrary diffusion coefficient in the very high energy regime, 
where the energy losses of particles dominate over their ac- 
celeration. In this regime, it is convenient to present the so- 
lution for the steady-state energy distribution of electrons 
in the following form: 



Ni.2 = A'i ! 2(x,p)exp[5 , i i2 (a:,p)] 



(3) 



where indices 1,2 correspond to the upstream and down- 
stream regions, respectively. In the asymptotic regime, 
when losses dominate over the acceleration, the function 
S(x,p) obtains large negative values. According to the stan- 
dard asymptotic method we assume that S(x,p) is propor- 
tional to b(x,p). After substitution of Eq. (3) into Eq. (1) 
and keeping all terms which are proportional to b 2 and 6, 
one finds 



Di, 
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dx 
d 2 S 1 
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dx 2 
dS h2 



2D 



dp 

dK h2 dS 1>2 
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(4) 



dx 



p 2 dp 



P 2 bi, 2 K 1:[ 



Denoting by 

S a (p) = Si(p)\ x=0 = S 2 (p)\ x=0 , 
K (p) = ^i(p)L =0 = K 2 {p)\ w=0 , 



(5) 

(6) 
(7) 
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and performing the same expansion for the boundary con- 
dition given by Eq.(2), we obtain 



uipdSo _ dS 2 
7 S dp dx 

mp dK dK 2 
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7s OP 



dx 



dSi 

dx 

dK x 



x=0 



dx 



(8) 



(9) 



The derivatives of functions Si and S 2 on x and the 
derivative dSo/dp can be found from Eqs.(4) and (8): 



dSi 
dx 

dS 2 
dx 

dS 



x=0 




Uip V i-> 2 



dp u\p- V 



(10) 

(11) 

(12) 



The solution of the last equation is 

r p dn' / , , \ 2 

S o = - 1 l 4 + y/b 2 {p')D 2 {p')) .(13) 

u Jq p \ / 

Taking the first derivative of Eq.(4) on x, and using Eqs 
(10) and (11) we find 
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After substitution of these functions into Eq. (5), and 
using Eq. (9) , we arrive at the following ordinary differential 
equation 



dK . 

P— An 



D 2P J 



^iP^Vh 



(16) 



dp u 2 VhD 1 ~+VhD 2 ~ 

the solution of which can be written in the following form 



K (p) oc p 1 / 2 exp J dp- 



V&i-Di + ^b 2 D 2 



(17) 



Eqs. (3), (13) and (17) determine the asymptotic form of the 
electron spectrum at the shock front Nq (p) for an arbitrary 
diffusion coefficient. In the case of an energy-independent 
diffusion, it coincides with the exact solution derived by 
Webb et al. ((19841) . 

These expressions can be significantly simplified in the 
case of the same energy-dependence of the diffusion coeffi- 
cient and the energy lose rate upstream and downstream: 
b 2 = b = D 2 = D = kD\. In this case 



(18) 



Ko(p) oc y/b/p , 

and Eq. (3) is reduced to 



Nn cx 



■ exp 



7s 



p dv' 
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.(19) 



The integrations of F\ = J° oQ dxN(x,p) and F 2 = 
J °° dxN(x,p) give the integrated electron spectra F\(jp) 
and F 2 {p) in the upstream and downstream regions, respec- 
tively. They are of great practical interest in many astro- 
physical situations. At large momenta they are determined 
by the ratios of the electron spectrum at the shock front 
to the absolute values of derivatives given by Eqs. (10) and 
(11): 



Fi(p) = 



uip 



^j s b(p) 



F 2 (p) 



uip 



1 + V^?7a&(p) 



N (p) 



No( P ) 



(20) 



(21) 



Eqs. (19)-(21) describe the form of the spectrum at large 
energies. At small energies the losses are negligible and the 
electron spectrum is power-law with an index 7 S . In order to 
describe the transition between this two extreme regimes, 
Eq.(l) has been solved numerically using an implicit finite- 
difference scheme. 

2.1. The case of Boh m diffusion 

Below we consider the most interesting case which as- 
sumes that the diffusion of electrons proceeds in the so- 
called Bohm diffusion regime and the energy losses of elec- 
trons are dominated by synchrotron radiation. We present 
the diffusion coefficient in the form D = r)cr g /3, where 
r g = pc/qB is the gyroradius of particles, B is the magnetic 
field strength, q and m are the electric charge and mass of 
the electron, respectively. The factor 77 > 1 allows deviation 
of the diffusion coefficient from its minimum value 77 = 1 
(the nominal Bohm diffusion). The synchrotron loss rate av- 
eraged over the pitch angles is b(p) = 4g 4 i? 2 p 2 /9rn 4 c 6 . Since 
the Bohm diffusion and synchrotron losses are determined 
by the magnetic field strength, the parameter £ = k~ 2 , 
where k is the ratio of the magnetic field upstream to the 
magnetic field downstream, k = B\jB 2 . Eq. (19) shows 
that in the high energy cut-off region the spectrum has a 
"Gaussian" type behavior, Nq cx exp(— p 2 jp\), where 



^ = (1 

mc 



r i/2\-l. 



2.86- 10 s 



7 s r ? 1 /2(i + K i/2) ^3000 km 



7sA /2w } £/27 



B 



100 /iG 



-1/2 



(22) 



Here the shock speed u\ and the magnetic field in the down- 
stream region B 2 = B are normalized to 3000 km s _1 , and 
100 fiG, which are quite typical for shells of young SNRs. 

At p = po, the characteristic lifetime of electrons due to 
the synchrotron cooling is 



= Po/b( Pa ) = 8.6 7s (1 + K 1 / 2 )t ? 1 / 2 x 



3000 km s 



B 



100 ^G 



-3/2 



yr 



(23) 



The synchrotron losses have a strong impact on formation 
of the energy spectrum of electrons at energies for which 
this time is small compared to the characteristic dynamical 
times of the source, e.g. the age of the accelerator. If this 
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condition is fulfilled, the electrons with such energies are 
concentrated in the shock vicinity, and therefore the plane 
shock approximation is well justified. 

In the case of standard shock acceleration in the Bohm 
diffusion regime with k = 1 and 7 S = 4, we combined the 
asymptotic analytical solutions, obtained at low and very 
high energies, with numerical calculations performed in the 
transition region. This allows a simple analytical presenta- 
tion of the electron spectrum at the shock over the entire 
energy range: 



N (p) cx p' 4 [l + 0.66(p/p ) 5/2 ] 9/5 exp( V/Po) 



(24) 



The spectrum of electrons given by Eq.(24) is shown in 
Fig. 2. At low and large energies Eq. (24) coincides with 
the exact asymptotic solutions: N (p) cx p~~ A and N (p) cx 
p 1 ! 2 exp(— p 2 /pq), respectively. In the cutoff region, the 
spectrum is proportional to the product of two terms - 
a power-law term with a positive slop (p 1 ^ 2 ) and an ex- 
ponential term, exp(— P 2 /Po)- While the first term can be 
interpreted as a pile-ujo, it in fact does not sow-up because 
the second (exponential) term effectively cancels this fea- 
ture. As a result, an almost perfect power-law (p~ A type) 
spectrum is formed up to p ~ po, with a super-exponential 
cutoff afterwards (see Fig. 2). 

It is easy to relate, using Eq. (1) and boundary condi- 
tion given by Eq. (2) the integrated spectra in downstream, 
upstream and the spectrum at the shock: 



Po 



F 2 + k^Fi = l—N (p), and I 



Isbpo ' 



(25) 



This relation is valid for all energies. It implies that for k = 
1 (i.e. equal magnetic field strengths upstream and down- 
stream) the integrated spectrum F = Fj + F 2 oc Nqp^ 1 . 

The case of different magnetic fields in the downstream 
and upstream regions, more specifically, for the scenario 
when the magnetic field downstream is stronger than mag- 
netic field upstream by a factor of n^ 1 = y/Tl, is of practi- 
cal interest, in particular for nonthermal emission of young 
supernova remnants. This case, which corresponds to the 
increase of the isotropic random B-field at the shock with 
the "standard" compression factor a = 4, can be described 
by the following analytical approximations for the electron 
distributions at the shock front: 

N (p) cx (p/p )- 4 [l + 0.523(p/p ) 9/4 ] 2 exp( -p 2 /p 2 ) , (26) 
and in the upstream region 

Fi(p) cx 0.70l(p/p )- 3 ll + 1.7(p/pof} 5/6 exp(-p 2 /p 2 ) .(27) 

For the downstream, the integrated energy spectrum of 
electrons can be found from Eq. (25): 



Po 



F 2 (p)^N {p)l^-0.09F 1 (p) 



P 



(28) 



Note that the spectrum integrated over the upstream 
region, at small energies is very flat, F\ cx pN (p) cx p~ 3 . 
This has a simple explanation related to the spatial scale 
of the electron distribution; in the upstream region it is 
proportional to the diffusion coefficient. 



1 The possible appearance of a pile-ups in the electron spectra 
within the so-called "box" model of shock acceleration has been 
noticed by Drury et al. (|1999p . 



The integrated spectra of electrons in the upstream and 
downstream regions given by Eqs. (26)- (28) as well as re- 
sults obtained numerically are shown in Fig. 3. 

The dependence of the spectra given by Eqs. (24), (26)- 
(28) on the factor rj is expressed through the momentum 
Po- 



3. Synchrotron radiation 

The synchrotron emissivity of electrons is determined as 



Zirmc z 



(29) 



where lo c is the characteristic frequency of synchrotron ra- 
diation ui c = 1.5 qBp 2 /in 3 c 3 . The function R describes syn- 
chrotron radiation of a single electron in magnetic field with 
chaotic directions. Crusius & Schlickeiser (|1986p derived an 
exact expression for R(lo/lo c ) in terms of Whittaker's func- 
tion. With an accuracy of several percent R(oo/u) c ) can be 
presented in a simple analytical form 



R(w/u c ) 



1.81 cxp(— lo/lo c ) 
V (lo c /lo) 2 / 3 + (3.62/tt); 



(30) 



These functions are compared in Fig. 4. 

The energy flux (J (to) = u)dN la d/dbo) of synchrotron ra- 
diation produced at the shock is determined by the integral 
of the emissivity along the line of sight I and over the solid 
angle f2 J(oj) = J dldfte. This means that it is determined 
by the integrated spectrum F(p). 

The integration of Eq. (29), using Eqs. (24) and 
(30), results in the energy flux of radiation J (to) = 
J dpa(p,uj)exp(-g(p,uj)), where g(p,w) = uj/uj c +p 2 /pl, 
a(p,u>) cx p 2 F(p)R(u>/uj c ) exp(g(p, u>)). At large frequen- 
cies, when the argument of the exponent is large this 
integral is mostly determined by the momentum in- 
terval around the momentum which minimizes the 
function g(jp,uS). Thus at large frequencies the en- 
ergy flux of radiation can be written as J{uj) = 
a(p*,u) exp(— g(p*, uj))yj2-K/g"(p*,uj). Here g" is the sec- 
ond derivative of g on p. The momentum p* is given by the 
expression: 



p* 
mc 



= 9 



,3„3„ ,2 



1/4 



1 



Vl + K 1 / 2 



= 8.3 x 10'x 



(Tiw/lkeV) 1 / 4 / m \ 1/2 / B 



7s V V /4 Vl + k 1 / 2 V 3000km A 



100 fiG 



-1/2 



(31) 



Note that at the fixed u>, the characteristic frequency of 
the synchrotron radiation to c corresponding to p» , is smaller 
than lu. 

In the regime of Bohm diffusion the energy flux of syn- 
chrotron radiation has an asymptotic form 



J cx w 3 / 8 exp(— ^co/wq) 
where 



81 



MTDIC 



UJ 



16 (l + K 1 / 2 ) 2 VlU 2 



(32) 



(33) 



V.N. Zirakashvili, and F.Aharonian: Energy spectra of electrons accelerated in supernova remnants 




0.01 - 



0.001 - 



1000 



Fig. 2. The electron distribution at the shock front (solid 
line) and the spatially integrated spectrum of synchrotron 
radiation (dashed line) produced by electrons accelerated 
at the shock with compression ratio a = 4 and equal up- 
stream and downstream magnetic fields (k = 1). The spec- 
trum of synchrotron radiation obtained using the S- function 
approximation is also shown (dotted line). 
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Fig. 3. The integrated upstream (Fi; dotted line) and 
downstream (F 2 ; dashed line) electron spectra given by Eqs. 
(27), (28), as well as the overall F = Fi + F 2 spectrum 
(solid line) for the shock with compression ratio <r = 4 and 
the ratio of the magnetic fields downstream and upstream 
fc _1 = vTT. The spectra F\ and F 2 obtained numerically 
are also shown (thin solid lines). 



or for the corresponding energy of synchrotron photons 

2.2 keV i ai \ 2 1 

7, 



huj, 







t ? (1 + k 1 /2)2 V 3000 km s 



\ 2 16 



(34) 



At small energies the integrated downstream syn- 
chrotron spectrum is proportional to uj ~. In order to 
describe analytically the synchrotron spectrum over the en- 
tire energy range, which should coincide with the exact 
asymptotic spectra at low and high energies, we use nu- 
merical calculations performed in the intermediate energy 
range. In particular, in the case of "/ s = 4 and k = 1, we 
suggest the following approximate formula 



J(u>) OC 



LO 



1 + 0.46 



wo 



11/4.8 



exp(-A/w/w ) , (35) 



which provides an accuracy better than 10 percent. 

The energy spectrum of electrons at the shock N (p) 
(see Eqs. (24)), is shown in Fig. 2. Because the cooling 
length of electrons downstream the shock is inversely pro- 
portional to the momentum, the spatially integrated elec- 
tron spectrum appears steeper, in particular, in the case of 
k = 1, F(p) = Fi{p) + F 2 (p) cx N (p)/p. Note that below 
the cutoff region F(p) oc p~ 5 . The corresponding spectrum 
of synchrotron radiation of these electrons is shown in Fig. 2. 
At small frequencies, J(u>) oc lu^ 1 . 

The synchrotron spectrum obtained using the so-called 
S- function approximation is also shown in Fig. 2. It is as- 
sumed, that the synchrotron emission of the electron is con- 
centrated at the frequency w = 0.29u c . This corresponds to 
the function R{lo/lo c ) oc 3(0.29 — lo/lu c ). Obviously this ap- 
proximation does not provide a satisfactory fit. 

In the case of j s = 4, n = 1/vTT we found the following 
spectra produced upstream and downstream the shock: 



Ji(w) oc 0.175 



1 + 1.27 



3/4' 



1/2 



J 2 {LO) OC 




11/4 




(37) 



exp 



3.32- 



Note that the spectrum of synchrotron radiation produced 
upstream the shock at small frequencies is very flat, J(ui) = 
const, which reflects the hard electron spectrum in that 
region (see Eq.(27) and Fig. 3). 

The spectra of synchrotron radiation in the upstream 
and downstream regions given by Eqs. (36) and (37) are 
shown in Fig. 5. The synchrotron flux in the downstream 
region is significantly enhanced because of the larger (com- 
pressed) magnetic field. The corresponding spectra ob- 
tained numerically are also shown. 

We should note that the exponential terms in Eqs. 
(35), (36) and (37) are relatively slow functions (propor- 
tional to exp[— (w/wo) 1 / 2 ])- Moreover, these exponential 
terms are multiplied to power-law terms (lu/ujq) s with pos- 
itive indices, s > 0, which effectively compensate the ex- 
ponential terms at frequencies lo <~ uj - For that reason 
the parameter luq only formally can be treated as a cut- 
off frequency. In reality, the breaks (cutoffs) in the syn- 
chrotron spectra appear at much higher frequencies, namely 
at lo > lOwo- 

The integrated downstream electron spectrum F 2 in 
Fig. 3 and the corresponding synchrotron spectrum J 2 in 
Fig. 5 may be applied to the supernova remnants with ages 
lager than the cooling time given by Eq.(23), i.e. when the 
electron spectrum in the cut-off region is determined by 
the energy losses of electrons. Electrons with energies below 
cut-off energy are accelerated with power-law spectrum and 
advected downstream. It was also assumed that these elec- 
trons lose all their energy downstream that is of course valid 
for energies larger than the energy of the electron which has 
the cooling time comparable with the age of the remnant. 
For smaller energies a three-dimensional geometry of the 
shock wave should be taken into account. Consequently, 
the integrated downstream electron spectrum and the cor- 
^3gy"esponding synchrotron radiation spectrum appear flatter 
compared to the ones shown in Figs. 3, 5. 
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Fig. 4. Comparison of the function R(uj/lo c ) tabulated by 
Crusius and Schlickeiser (JT986) (dotted line) and the func- 
tion given by the approximate expression (30) (solid line) . 
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Fig. 5. The energy spectra of synchrotron radiation given 
by Eqs. (36), (37) and produced by electrons in downstream 
(dashed line) and upstream (dotted line) regions of the 
shock with compression ratio a = 4 and the ratio of the 
magnetic fields downstream and upstream = VdT. The 
corresponding spectra obtained numerically are also shown 
(solid lines). 

3.1. The radiation zones 

The asymptotic solutions derived in Sec. 2 can be used for 
study of coordinate-dependence of electron distributions. 
Near the shock front, the spatial distribution of electrons 
has an exponential dependence. The inverse exponential 
lengths in the upstream and downstream regions are given 
by Eqs. (10) and (11). Since the synchrotron radiation is 
produced mainly by electrons with momentum , the ex- 
ponential lengths of emissivity (29) in the upstream l u and 
downstream Id regions of the shock can be found from Eqs. 
(10), (11) and (31): 



k 
and 



9 / m 9 c 17 u 2 ri 

V3(l + kV2) = W 2 £ 6 7? 



-5/4 



(38) 



(39) 



In particular, for 7 S = 4 we have 
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(40) 



for the radiation width in downstream, and a factor of /t~ 5 / 4 
larger in upstream. Note that the above equations are cor- 
rect only for large frequencies, oj > ojq. 

The radiation width in downstream (40) is relatively 
small. This permits to explain the narrow X-ray fila- 
ments observed in SNRs and to estimate the magnetic field 
strength (see e.g. Berezhko et al. 120051 Vdlk et al. [5005)) . 

4. Inverse Compton scattering 

In the isotropic photon field, the emissivity of the inverse 
Compton radiation of an electron with Lorentz factor 7 is 
described as (Blumenthal & Gould [T970) 



eio(w) = 2hcr 2 
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0J C 



1 + 2 In ■ 
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Here uj c — Auj lr ) 2 (l 



fhUJ 



r e = q 2 /mc 2 is the classi- 
cal electron radius and n(u>) is the energy distribution of 
the target photons. Hereafter we assume that the target 
photons are described by a Planckian distribution 



n(u>) 



W 2 /(C 3 7T 2 ) 



exp(ftw / kbT) — 1 



(42) 



where T is the temperature of radiation and Kb is the 
Boltzman constant. It is easy to show, using Eqs. (41) 
and (42), that at large frequencies oj >> ^kbT/H the IC 
emissivity in the Thompson regime behaves as eic(w) oc 
oj exp(— huj / Aj 2 kbT) . At very high energies the integra- 
tion over the electron energy distribution at the shock and 
the use of the method described in the Section 3, give 
the following asymptotic form of the IC spectrum in the 
Thompson limit: 



J IC cx w 7/8 exp(- V^M) > 
with 

27 KBTu\m 2 c 2 

LU b - 



(43) 



(44) 



2?7(l + K i/2) 2 hj 2 q 3 B 
or, in terms of energy of gamma-rays, e b = hu b , we have: 
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(45) 



As before, it is assumed that synchrotron losses dominate 
over IC losses, i.e. the energy density of the background 
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Fig. 6. The spectrum of IC radiation given by Eqs. 
(47), (48) and produced by electrons accelerated down- 
stream (dashed line) and upstream (dotted line) the shock 
with compression ratio a — 4 and the ratio of the mag- 
netic strength downstream and upstream \f\\. The sum 



J lc = Jl c + Jl c is shown by the solid line. It is as- 
sumed that IC scattering proceeds in the Thompson limit, 
&kn = 0, and that synchrotron losses dominate over the IC 
loses. The corresponding spectra obtained numerically are 
also shown (thin solid lines). 



radiation C/ rad << B 2 /8tt. In the case of 2.7 K CMBR this 
implies B > 3/iG. 

Using the exact asymptotic spectra at low and high en- 
ergies, and the results of numerical calculations at interme- 
diate energies, we can describe the broad-band energy flux 
of the inverse Compton (Thompson) scattering by simple 
analytical expressions. In particular, for k = 1 one has 



tIC / \ w & 
J (uj) oc — 

to 



0.36 | — 

,Ub 



0.7 



15/5.6 




For the "standard" scenario of shock acceleration with 
7 S = 4 and k = 1/ VTT, we obtain the following IC spectra 
produced upstream and downstream the shock: 



J\ C {uj) oc 0.20 



D.7.". i — y /s 




exp [-J- , (47) 
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The spectra given by Eqs. (47) and (48) are shown in Fig. 6. 
They are derived under the assumption that the Compton 
scattering occurs in the Thompson regime. However, in 
many cases one should take into account the Klein-Nishina 
effect. The importance of this effect in the case of black- 
body radiation field is determined by the parameter 
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Fig. 7. The spectra of IC radiation produced by electrons 
accelerated at the shock with compression ratio a — 4 and 
the ratio of the magnetic strength downstream and up- 
stream "v/ll for different values of bxN- bxN — (solid 
line), bxN — 0.1 (dashed line), bxN — 0.3 (dotted line), 
bKN = 1.0 (thin solid line), bKN = 3.0 (thing dashed line) 
and bKN = 10.0 (thin dotted line). The case of dominating 
synchrotron losses is considered. 



The Thompson limit corresponds to b^N << 1> while the 
Klein-Nishina effect becomes significant already at bj<N > 
0.1. The IC spectra of electrons calculated for different val- 
ues of the parameter bKN are shown in Fig. 7. 

As long as the energy losses of electrons are dominated 
by synchrotron radiation, the Klein-Nishina effect makes 
steeper the spectrum of IC gamma-rays. In particular, it 
moves the region of the spectral cutoff to lower energies. 
Below the cutoff region, the energy spectrum of IC gamma- 
rays still can be described by a (quasi) power-law. The 
Klein-Nishina effect changes the slope of the gamma-ray 
spectrum significantly, even for relatively small values of 
the parameter bKN- 



5. Discussion 

The detection of nonthermal X-ray emission of synchrotron 
origin clearly demonstrates the existence of multi-TeV elec- 
trons in different astrophysical source populations like shell 
type SNRs - pulsar wind nebulae, microquasars, small 
and large scale jets in AGN, etc. Given the severe syn- 
chrotron losses of electrons, which increase with energy as 
dE/dt oc E 2 , the acceleration efficiency in these objects 
should be very high in order to boost the electrons to en- 
ergies well beyond 1 TeV. In particular, in young SNRs 
with typical shock speeds u\ ~ 2000 — 3000 km/s, which 
are established as prominant nonthermal X-ray emitters, 
the diffusive shock acceleration of electrons should proceed 
in the regime close to the Bohm diffusion limit, otherwise 
the synchrotron spectrum would break before achieving the 
X-ray domain. 

Although the basic concepts of diffusive shock accel- 
eration are deeply studied and well understood (see e.g. 
Malkov and Drurv [2"00ip . many important details remain 
unexplored. This concerns, in particular the radiation fea- 
tures of the shock accelerated electrons. Over the recent 
years, several numerical calculations of high energy radia- 
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tion by particles accelerated in young SNRs have been per- 
formed with emphasis on the non-linear effects caused by 
protons and nuclei on the shock structure (see e.g. Ellison 
et al. 120001 Berezhko and Volk l2"005|l . 

The shock accelerated electron spectra have been stud- 
ied analytically (e.g. by Webb et al. I1984p , but only for the 
case of energy-independent diffusion coefficient. Therefore 
the results of these early works generally cannot be ap- 
plied to the most interesting sources with synchrotron X- 
ray emission. These sources, as noticed above, require dif- 
fusive shock acceleration of electrons in the regime close to 
the Bohm diffusion which is characterized by a diffusion 
coefficient D(E) oc E. 

The aim of this paper was a fully analytical treatment 
of the energy spectra of shock-accelerated electrons and 
their synchrotron and inverse Compton radiation for an 
arbitrary energy-dependent diffusion coefficient. In partic- 
ular, under the assumption of dominance of synchrotron 
losses, we derived exact asymptotic solutions in the high 
energy region. For the energy-dependent diffusion coeffi- 
cient written in a rather general form D oc p a , the elec- 
tron spectrum contains an exponential term exp[— (p/po)^] 
with f3 = a + 1 (see Eq.(19)). Thus, only in the case of 
energy-independent diffusion (a = 0) the electron spec- 
trum is characterized by an exponential cutoff. In the most 
interesting case of Bohm diffusion, the spectral cut-off is 
of Gaussian type, exp[— (p/po) 2 ]- The corresponding cut- 
off in the spectrum of syncrotron radiation appears much 
smoother (see also Fritz 1989]for discussion of this effect). It 
is described by an exponential term exp[— (w/wo)^' 2 ^'] = 
exp[— (lo / uJo)( a+1 ^ ( a+ ^\. For the energy independent diffu- 
sion coefficient we have J{lo) oc exp[— (lo/luq) 1 ^ 3 ], while in 
the case of Bohm diffusion J(w) oc exp[— (cj/ujq) 1 ^ 2 ]. 

We should note that the exp[— (lh/ujq) 1 / 2 } type behav- 
ior of the synchrotron spectrum in the cutoff region for- 
mally is similar to the one derived for the generally as- 
sumed "power-law with exponential cutoff' electron spec- 
trum, No(p) oc p~ r exp(— p/po), and using the 5-function 
approximation for calculation of the synchrotron radiation 
(see e.g. Reynolds 11998ft . However, this is simply a coinci- 
dence. First of all, the ^-functional approximation leads to a 
wrong synchrotron spectrum; the exponential cut-off in the 
electron spectrum results in the synchrotron spectrum with 
a exp(— w 1 / 3 ) type cut-off. More importantly, the spectra of 
shock accelerated electrons given by Eq.(26)-(28) are quite 
different from the simplified "power-law with exponential 
cutoff" assumption. Correspondingly, because of the depen- 
dence u> c oc p 2 B, the effect on the spectrum of synchrotron 
radiation appears even more pronounced (see Fig. 2). 

Exact solutions for electrons are possible only at low 
and high energies, namely in the energy intervals where (i) 
the radiative losses of electrons are negligible, and (ii) the 
radiative energy loss rate exceeds the acceleration rate. The 
indices of the power-law factors in the asymptotic solutions 
(19), (32) and (43) significantly differ from the power-law 
indices of the electrons, as well as the synchrotron and IC 
spectra below the cut-off regions. In order to connect these 
solutions we performed numerical calculations in the inter- 
mediate energy intervals. This allowed us to obtain, within 
accuracy better than 10 percent, simple analytical presen- 
tations for the spectra of electrons and their synchrotron 
and inverse Compton radiation components. 



Eqs. (35)-(37) and Eqs. (46)-(48) properly describe the 
relative contributions of radiation from the regions up- 
stream and downstream the shock. At the same time, they 
do not allow direct calculations of the absolute synchrotron 
and inverse Compton luminosities. However, the normali- 
sation constants can be easily calculated from the total en- 
ergy release through the synchrotron and inverse Compton 
channels. 

It should be noted that in this paper the momentum po 
in the electron spectrum, and the corresponding frequencies 
too and cob which characterize the positions of cut-offs in the 
spectra of synchrotron and IC radiation components, are 
found from the exact asymptotic solutions of kinetic equa- 
tions. Interestingly, the results appear quite close, within 
a factor of 2, to the estimates derived from the condition 
"acceleration rate = synchrotron cooling rate". However, 
as it follows from Fig. 2-6, the real cut-offs appear, in fact, 
at much higher frequencies. In the case of synchrotron ra- 
diation, the spectrum continues without an indication of 
a break or a cutoff up to 20ojq. The spectrum of inverse 
Compton gamma-rays turns down also at u> ~ 20wf, (in 
the case of constant magnetic field), or it continues up to 
~ 50wb (in the case of different magnetic fields upstream 
and downstream) . The more effective extension of the hard 
spectrum of IC radiation is explained by the larger spatial 
scale of electron distribution upstream the shock. This is 
not important for the synchrotron radiation the contribu- 
tion of which from the upstream region is suppressed by an 
order of magnitude. On the other hand, because of the ho- 
mogeneous distribution of the target radiation field, 2.7 K 
CMBR, the IC gamma-ray contribution at highest energies 
from the upstream region dominates over the contribution 
from the downstream region. 

The results obtained in this paper can be used to fit the 
X-ray spectra observed from young supernova remnants. 
The value of the parameter loq found from the fitting pro- 
cedure should be compared with the theoretical value given 
by Eq.(33). For the known shock velocity of a given super- 
nova remnant, this would allow an important estimate of 
deviation of the diffusion regime from the extreme Bohm 
limit. The dependence of the energy spectra on the factor 
•q is expressed via parameters p^ 1 u)q, lo\,. If the width of 
X-ray filaments is known, the magnetic field strength may 
be found from Eq. (40). All this in turn would provide reli- 
able estimates of the maximum energy of protons and nuclei 
accelerated in supernova remnants. 

Finally, we note that the results of this paper are ap- 
plicable for nonrelativistic shocks with velocities up to 
Mi ~ (0.1 — 0.3)c. The derived analytical approximations 
do not take into account possible nonlinear effects which 
can somewhat modify the spectra of electrons at low and 
high energies. 
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